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We investigate the singlet-triplet relaxation due to the spin-orbit coupling together with the 
electron-phonon scattering in two-electron multivalley silicon single quantum dots, using the exact 
diagonalization method and the Fermi golden rule. The electron-electron Coulomb interaction, 
which is crucial in the electronic structure, is explicitly included. The multivalley effect induced by 
the interface scattering is also taken into account. We first study the configuration with a magnetic 
field in the Voigt configuration and identify the relaxation channel of the experimental data by Xiao 
et al. [Phys. Rev. Lett. 104, 096801 (2010)]. Good agreement with the experiment is obtained. 
Moreover, we predict a peak in the magnetic-field dependence of the singlet-triplet relaxation rate 
induced by the anticrossing of the singlet and triplet states. We then work on the system with 
a magnetic field in the Faraday configuration, where the different values of the valley splitting 
are discussed. In the case of large valley splitting, we find the transition rates can be effectively 
manipulated by varying the external magnetic field and the dot size. The intriguing features of the 
singlet-triplet relaxation in the vicinity of the anticrossing point are analyzed. In the case of small 
valley splitting, we find that the transition rates are much smaller than those in the case of large 
valley splitting, resulting from the different configurations of the triplet states. 

PACS numbers: 73. 21. La, 71.70.Ej, 72.10.Di 



I. INTRODUCTION 



Spin-based qubits in semiconductor quantum dots 
(QDs) are believed to be a promising candidate for 
scalable quantum computation.^ Among different kinds 
of QDs, GaAs ones have been extensively investigated 
in the past decadei^"— As reported, the spin deco- 
herence, which is essential to know for genuine appli- 
cations in such systems, is limited by the hyperfine 
interactioni^"— and the spin-orbit coupling (SOC) ^^i^° 
together with the scatteringi^"— Recently, much at- 
tention has been given to silicon due to its outstanding 
spin coherence properties i^"— Specifically, the hyperfine 
coupling strength in natural silicon is two orders of mag- 
nitude weaker than that in GaAs^S and can be further 
reduced by isotopic purification. jjj addition, the Dres- 
selhaus SOC^^ is absent in bulk silicon because of the 
existence of the bulk inversion symmetry. Although the 
interfaces of a confined system can introduce an interface 
inversion asymmetry (IIA),^^"'*'^ the SOC due to this ef- 
fect is still very small. Moreover, the absence of the piezo- 
electric interaction makes the electron-phonon scattering 
much weaker compared to that in III-V semiconductor 
QDs.'^^ All these features suggest a long spin decoherence 
time in silicon QDs, which is of great help in realizing the 
operation of logic gates and the storage of information. 
Furthermore, as an indirect gap semiconductor, silicon 
has sixfold degenerate minima of the conduction band. 
This degeneracy can be splitted by strain or confinement 
in quantum wells into two parts: a double-degenerate 
subspace of lower energy and a fourfold-degenerate sub- 
space of higher energy. The presence of the interfaces can 
further lift the twofold degeneracy by a valley splitting 



energy, which is strongly dependent on the size of the 
confinement structure4^i^ Moreover, the correlation ef- 
fects in silicon are much stronger than those in GaAs due 
to the enhanced electron-electron Coulomb interaction, 
thanks to the reduced kinetic energy due to the larger 
effective mass. Thus, the physics in silicon is expected to 
be richer. 

Very recently, spin-qubits utilizing the singlet- 
triplet (ST) states in silicon QDs have been actively 
investigated^ir— Culcer et al^ analyzed the feasibility 
of initialization and manipulation of ST qubits in dou- 
ble QDs, concentrating on the multivalley effect. With a 
large valley splitting, the exchange coupling was explic- 
itly investigated by Li et al.M- However, to the best of our 
knowledge, study on the ST relaxation in silicon QDs is 
rather limited,^'^i2&'^ Prada et alP calculated the ST 
relaxation using the perturbation method with the low- 
est few levels, which has been shown to be inadequate to 
study the ST relaxation time?^ Moreover, the Coulomb 
interaction was not explicitly calculated in their work. 
However, the strong Coulomb interaction together with 
the SOC are of critical importance to determine the ener- 
gies and wave functions of the eigenstates in QDs. There- 
fore, the diagonalization approach with a large number 
basis functions is required to guarantee the convergence 
of the energy spectrum and the ST relaxation rates 
In the present work, we calculate the ST relaxation in 
silicon single QDs by explicitly including the Coulomb 
interaction and the multivalley band structure as well 
as the SOCf^i^ which is the key of the ST relaxation 
mechanism discussed in this work. In the calculation, we 
employ the exact diagonalization method and calculate 
the ST relaxation rate from the Fermi golden rule^i^ We 
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first calculate the ST relaxation rate in silicon QDs with 
a parallel magnetic field (i.e., the Voigt configuration). 
Our theory successfully explains the recent experiment, 
by Xiao et alM- and suggests that the measurement cor- 
responds the relaxation of the lowest singlet, with the 
dominant channel being the one associated with the low- 
est triplet. We further predict a peak in the magnetic- 
field dependence of the ST relaxation rate, resulting from 
the large spin mixing at the anticrossing point between 
the singlet and triplet states. Then we investigate the 
perpendicular magnetic-field (the Faraday configuration) 
dependence of the ST relaxation rate with different val- 
ues of the valley splitting. In the situation of large valley 
splitting, the lowest singlet and three triplet states are 
all constructed by the lowest valley state. We find that 
the transition rates can be effectively manipulated by the 
magnetic field and dot size. We also find intriguing fea- 
tures in the vicinity of the anticrossing points. Moreover, 
we compare the relative contributions of the intravalley 
transverse acoustic (TA) and longitudinal acoustic (LA) 
phonons to the transition rates. In the case of small 
valley splitting, the eigenstates of the lowest two valleys 
contribute. We find the ST relaxation in this case is much 
slower than that in the large valley splitting one. 

This paper is organized as follows. We set up the model 
and give the formalism in Sec. II. Then in Sec. Ill, we 
utilize the exact diagonalization method to obtain the 
energy spectrum and calculate the ST relaxation rates. 
Both parallel and perpendicular magnetic-field depen- 
dences of the ST relaxation rates are studied. The behav- 
ior of the transition rates in the vicinity of the anticross- 
ing points is also discussed in this section. We summarize 
in Sec. IV. 



II. MODEL AND FORMALISM 

In our model, we choose the lateral confinement as 
Vc{x^y) = imtu;o(x^ -|- y^), with rrit and cjq represent- 
ing the in-plane effective mass and the confining poten- 
tial frequency.—"^ The effecti ve diameter can then be 
expressed as do = ^hn /mtUJo- In the growth direction 
[001], Vz{z) is applied within the infinite-depth well po- 
tential approximation. Then the single-electron Hamil- 
tonian with magnetic field B = B±z -\- Byx is described 
by 

ge = P-^ + Py\ E^ + v{v) + (P) + Hz+H,, (1) 

with rUz representing the effective mass along the z- 
direction. V{r) = Vc + Vz and P = -ihV + (e/c)A, 
with A = {-yB±,xB±,2yBn)/2. H^o describes the SOC 
Hamiltonian, including the Rashba term-° due to the 
structure inversion asymmetry and the term due to the 
IIA4ii"i^ Then, one obtains 



where oq and bo stand for the strengths of the Rashba 
and IIA terms, respectively. The Zeeman splitting is 
given by Hz = ^gi^B{B±(7z + B^^a^) with g being the 
Lande factor. Since the four in-plane valleys are sepa- 
rated from the two out-of-plane ones by a large energy 
gain, only the two out-of-plane valleys are relevant in the 
calculation. Hv in Eq. ([T]) describes the coupling^^^ be- 
tween the valleys lying at ±(A:o) along the z-axis with 
(ko) — 0.85(27r/asi). Here, asi = 5.43 A is the lattice 
constant of silicon. Throughout this paper, one uses the 
subscripts "z" and "z" to denote the valley at (fco) and 
the one at —(fco), respectively. 

By solving the Schrodinger equation of the Hamilto- 
nian Hq = ^{Px^ + Py^) + 2^Pz'^ + ^(r), one obtains 



the eigenvalues^!^ 

Er,ln, = Kl{2n + \l\ + 1) + hlLOB + 



Sniza'^ 



(3) 



where fl — ^/luq'^ + lob'^ and lob — eB^/(2mt). a repre- 
sents the half-well width. Here, n — 0, 1, 2, ... is the ra- 
dial quantum number and / = 0,±1,±2,... represents the 
azimuthal angular momentum quantum number. The 
index denotes the subbands resulting from the con- 
finement along the growth direction. The corresponding 
eigenfunctions read 



^sin[I^(z-Ka)], |z|<a 



0, 



otherwise. 



(4) 



with Nr.,1 = {a'^n\/[K{n+ 1^1)!]}^/^ and a = ^mt^jh. 

l}n is the generalized Laguerre polynomial. The wave 
functions in different valleys can then be expressed as 



Fr,. 



F(r), with M2,j(r) representing 
the lattice-periodic Bloch functionsi^ Here, we neglect 
the orbital effect of the parallel magnetic field by consid- 
ering the strong confinement along the growth direction. 

One can demonstrate that the overlap between 
the wave functions in different valleys is negligi- 
bly small, therefore only is considered to con- 
tribute to the intervalley coupling in the present 
work. However, there still remain some controver- 
sies over the valley coupling nowadays.^^^'^^ In this 
work, we take (C/nJ-^vlC'Z'n,) = A^^(5„„'(5;z' and 
(CinJ-f^vlC'/'nJ = ^nJnn'Sw, accordiug to Ref.lH. 
Here, only the coupling element between the states with 
identical is given, since only the first subband is in- 
cluded in our calculation while the others are neglected 
due to the much higher energy. Including this interval- 
ley coupling, the eigenstates become = ^((^L„ ± 



with eigenvalues 
these equations. 



En 



±|Ai I. In 



Hso = aoiPxCTy - Py<7x) + bo{-Px<7x + PyCTy 



(2) 



A 



Arnztt^ ' 
VyHz'^TT^h'^ cos(2fcoa) 
AvrLzO? 



(5) 
(6) 



3 



with Vv standing for the ratio of the vahey couphng 
strength to the depth of quantum well4^ 

For a two-electron QD, the total Hamiltonian is given 
by 

Htot = (^^c + + He) + Hp + H^p + Hip. (7) 

Here, the two electrons are labeled by "1" and "2". 
The electron-electron Coulomb interaction is described 



by He 



with K, representing the relative 



4ireo«;|ri-r2| 

static dielectric constant. Hp = X^qA S^^qAflq^aqA rep- 
resents the phonon Hamiltonian with A and q denot- 
ing the phonon mode and the momentum respectively. 
The electron-phonon interaction Hamiltonian is given by 

We construct two-electron basis functions in the forms 
of either singlet or triplet based on the the single-electron 
eigenstates. For example, we use two single-electron spa- 
tial wave functions \nilinzinyi) and |n2/2'^z2?^D2) (de- 
noted as |iVi) and |A^2) for short; n.„ — ±) to obtain the 
singlet functions 

(8) 



and the triplet functions for Ni ^ Ni 



IT 



\T, 



(H)\ 



1 

1 



(|A^iA^2)-|7V2iVi))®|tt), 



(9) 



(i7Vi7V2) - |iV27Vi)) ® (in> + 1 ;t»(io) 



1 



(11) 



Here, the spatial wave functions of the first and second 
electrons in \NN') are denoted as N and N' in sequence. 
The superscript (S) denotes the valley configuration of 
each state. We define S = ± for the valley indexes of 
single electron states n„i — n^i — ±, and 5 = m for 



Then, one can calculate the matrix elements of the 
Coulomb interaction, which can be expressed by 



{N,N2\Hc\Kn;) ^ ^^^<S,^+,^,,+,, (12) 



where the superscripts 7^ and 7- run over the two val- 
leys, z and z, with yy^. = 1 and rfj^ = —t]^ — 1. G 
is given in detail in Appendix. One also calculates the 
SOC and Zeeman splitting terms, hence obtains the two- 
electron Hamiltonian, i.e., the terms in the bracket in 
Eq. ([7]). Then, the two-electron eigenvalues and eigen- 
functions can be obtained by exactly diagonalizing the 
two-electron Hamiltonian. We identify a two-electron 
eigenstate as singlet (triplet) if its amplitude of the sin- 
glet (triplet) components is larger than 50 %. 

From the Fermi golden rule, one can calculate the tran- 
sition rate from the state |i) to |/), due to the electron- 
phonon scattering, 

^^^f = YEl^^q>l'l(/lxN)n^qA'5(e/-e»-?l^qA) 

qA 

+ {nqx + l)6{ef-e, + nujqx)], (13) 

in which x(q, ri, r2) = e''^ ""^ -I- e^'^ '"^ and Uqx stands for 
the Bose distribution of phonons. In our calculation, the 
temperature is fixed at K. Thus n^x = and only the 
second term contributes. 



III. NUMERICAL RESULTS 

Since the piezoelectric interaction is absent in silicon — 
and the energy difference between the initial and fi- 
nal states discussed here is much smaller than the en- 
ergies of the intervalley acoustic phonon and the opti- 
cal phonon,— one only needs to calculate the intravalley 
electron-acoustic phonon scattering due to the deforma- 
tion potential. In the present work, both the TA and 
LA phonons are included. The corresponding matrix 
elements are M|ij^jj.3^q = hDjQ^ / {2dnfj^intra,Q) with 
/3=LA/TA standing for the LA/TA phonon mode. Here, 
we take the mass density of silicon d = 2.33 g/cm'^<^ The 
deformation potentials for the LA and TA phonons are 
Dla = 6.39 eV and Z^ta = 3.01 eV, respectively.^o The 
phonon energy r2^jntia,Q = VjsQ with sound velocities 
WLA = 9.01 X lO^cm/s and vta = 5.23 x lO^cm/s^ In 
our calculation, we take rrit = 0.19mo and = 0.98mo 
with mo being the free electron massi^ The Lande factor 
g = 2^ the ratio K = 7.2 x lO'^m.i^ In the calcula- 
tion, we employ the exact diagonalization method with 
the lowest 1516 singlet and 4452 triplet basis functions to 
guarantee the convergence of the energies and the tran- 
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sition rates. 

One finds that the eigenstates composed by the two- 
electron basis functions with single valley state "— " are 
almost independent of those constructed by the ones 
with single valley state "+" and two valley states "— " 
and "+" . On the one hand, there is nearly no cou- 
pling between them due to the negligibly small interval- 
ley Coulomb interaction'^^ and overlap between the wave 
functions in different valleys. One can also demonstrate 
that the elements of the SOCs between the states with 
different valley indices vanish when only the first subband 
is included, regardless of the coupling strengths ag and 
60. On the other hand, the transition between them is 
almost forbidden because (/|xl*) in Eq. ((T3)) is strongly 
suppressed thanks to the large intervalley wave vector 
(2fco) from the difference of the phases between different 
valleys. Therefore, the eigenstates are divided into three 
independent sets based on the valley indexes. It is noted 
that the energy of the eigenstate with valley configuration 
"— " is smaller than the corresponding levels with valley 
configurations and "m" due to the contribution of 
the valley splitting. 



A. PARALLEL MAGNETIC-FIELD 
DEPENDENCE 

Very recently, Xiao et ali^ measured the ST relaxation 
time in Si/Si02 QDs under a magnetic field parallel to 
the interface of the heterostructure. They reported that 
the ST relaxation time only slightly fluctuates around 
5 ms when the magnetic field increases from 2 to 4.5 T. 
In the experiment, the orbital level spacing is observed 
to be ~ 0.4 meV, corresponding to the effective diameter 
of the QD do = 56 nm. However, some parameters such 
as the effective well width, the strengths of SOCs and 
the valley splitting are unavailable. Moreover, the chan- 
nel of the relaxation process is not identified because of 
the uncertainty of the exact excited states spectrum in 
the experimenti^ Here, we take advantage of our model 
to clarify the experiment situation. In the calculation, 
we assume the magnetic field along x-direction and take 
the relative static dielectric constant k = 7.942. Since the 
valley splitting is strongly dependent on the effective well 
width according to Eq. (6), it is difficult to determine the 
energy spectrum without the knowledge of the exact well 
width. For a large valley splitting, the lowest levels are 
all constructed by the states with the single valley in- 
dex "— " , and the energy difference between the adjacent 
levels is determined solely by the orbital level spacing 
and Zeeman splitting approximately. Therefore, the re- 
laxation rate of each excited state can be calculated to 
identify the relaxation channel in the experiment. How- 
ever, the lowest levels become more complicated for a 
small valley splitting because new levels with the valley 
index "+" become relevant. Fortunately, as said above, 
the inclusion of the states with valley configuration 
or "m" has no observable influence on the relaxation of 
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FIG. 1: (Color online) (a) The lowest few energy levels vs. 
parallel magnetic field By in a single QD. Note that each 

of \T^~^) and prl"^) {i = 1,2,3,6) is double-degenerate, (b) 
Relaxation rates vs. the magnetic field. The red dots stand for 
the experimental data. S'~' "p*^"' (i = 1-3) denotes the sum 

of the relaxation rates from 15*^"') to the two degenerate 'T^^ 
levels. In the calculation, 2a = 4.344 nm and do = 56 nm. 



the states with the valley configuration "— " . In the fol- 
lowing, we first study the large valley splitting case. We 
take 32 monoatomic layers of silicon along the growth 
direction, corresponding to the well width 2a = 4.344 nm 
(2|Ai J = 0.83 meV). The strengths of the Rashba SOC 
and HA term are used as fitting parameters. We first cal- 
culate the energy spectrum since it is weakly dependent 
on the strengths of the SOCs. The lowest few levels, de- 
noted as IS'(-)), |T(j"^) and \'t[~^) (spin down) {i = 1-6) 
according to their major components, are plotted as func- 
tion of the magnetic field in Fig.fTfa). As the magnetic 
field increases, the energies of l^^^^) and \T^~^) keep in- 
variant while those of {i — 1-6) decrease due to 
the Zeeman splitting. The major component of 15*^^') is 
constructed by the single-electron states |01)|0 — 1) ac- 
cording to Eq. ([5]), and those of the triplet states |*t1 •*) 
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and l^rr') are given by |01)|0- 1) and |00)|10) follow- 
ing Eq. (dH), respectively. We find \'t[~^) {i = 1,2,3,6) 
and \Tq~^) are all double-degenerate. The major com- 
ponents of the two levels of |^t1"^), \^t[^^), I^t't^), 

and I^tI ■*) are composed by the single-electron states 
|00)|0±1), |0±1)|0±2), |00)|0±2) and |10)|1±1) in 

sequence. |Tq ') is mainly constructed by the basis func- 
tion involving |00)|0 ± 1) also. Here, we only retain the 
quantum numbers n and I for short, because other quan- 
tum numbers of these single-electron states are all the 
same. 

We then calculate the relaxation rates of these states 
due to phonon emission. Due to the low temperature in 
the experiment,"^** the relaxation rate at zero temperature 
can well represent the experimental data. We find that if 
one takes the Rashba SOC strength qq = — 2.09m/s and 
the IIA term strength bo — — 10.44 m/s, the total relax- 
ation rate of the state IS**-"^) fits the experimental data 
pretty well as shown in Fig.[TJb) (from 2 to 4.5 Tesla). 
The relaxation rates of other levels can not recover the 
experiment results. Specifically, the relaxation rate of 
\Tq presents a peak at ~ 3.45 T (not shown in the 

figure) and those of \^T^^) (i = 2-6) relax too fast (in the 
magnitude of ^ Ins). Therefore, we conclude that the 
experimental data might correspond to the lifetime of the 
smg let |S'(")). The rates of the major relaxation channels 
of IS'(-)) (involving S^") ^tL"\ i = 1-3) are also plot- 
ted in Fig.[ljb). Interestingly, the calculation predicts a 
peak of the total relaxation rate at By ~ 6.5 T, which 
should be checked by future experiments. Moreover, one 
also finds the significant increase of the total relaxation 
rate by increasing the magnetic field in the small mag- 
netic field regime, i.e., below 2 T. Such rich magnetic- 
field dependences can be understood as follows. From 
the figure, we find that the dominant relaxation channel 
is the one from |S'(-)) to I^tI ^). In the small magnetic 
field regime, the energy of the phonon emmision of this 
channel (corresponding to the energy difference between 

l^^")) and I^tI"-*)) is small and linearly increases with 
the magnetic field, which lead to the significant enhance- 
ment of the transitioui^iii However, the transition rate 
becomes insensitive to the phonon energy since the value 
of (/IxN) in Eq. is suppressed for a large phonon 
momentum, then the transition rate only slightly varies 
beyond 2 T. This picture can be also used to understand 
the feature of the relaxation rates between IS'*-"-') and 
{i = 2,3) far away from the peak. The peak at 
B\\ ^ 6.5 T, where the triplet state |^rl ■*) intersects the 
singlet state jS'^"^), results from the strong coupling be- 
tween them due to the SOCs. To ease further discussion, 
one denotes the total angular momentum and spin states 
as L = h+h and {S,Sx), respectively, with Sx represent- 
ing the X component of the total spin S. By neglecting 
the terms with ax in Eq. ^ due to its smaller magni- 
tude compared with the Zeeman splitting, one obtains 



the SOC Hamiltonian 



F,o = [ao{P+ + P-) - ibo{P+ - P-)] {S+ + S-)/K 

(14) 

with = Sy ± iSz- The ladder operations and 
change L and Sx by one unit, respectively. There- 
fore, a state with (L, Sx) can couple with the one with 
(L ±l,Sx± 1) or (L ±l,SxT 1) for both the Rashba 
and IIA terms. From the major components of the two- 
electron eigenstates, the quantum numbers {L,Sx) of 
I^tI"') and jS'^")) are (±1,-1) and (0,0), respectively. 

It is obvious that \^T^^) directly couples with IS"^"-*) 
through the SOCs. As a result, there is an energy gap 
(too tiny to pick up in the figure) at the intersecting 

point between |*^t1"^) and |S'(~'), which means an anti- 
crossing event occurs. In the vicinity of this anticross- 
ing point, the wave function of IS"^"-*) contains a large 
amount of the spin-down triplet component, which en- 
hances the spin relaxation process. One notices that the 
intersecting point between pTl and IS*'^') is also an 
anticrossing point. However, the coupling between these 
states is indirect and small, hence only slightly affects the 
ST relaxation. Other intersecting points between |S'*-~-*) 
and prl {i — 3-5) are just simply crossing points. 
Moreover, one finds that the peak of the relaxation rate 
of |Tq ') at B|| ~ 3.45 T refiects the anticrossing behav- 
ior between IT^"^) and \^t[^^) {i = 3-5) and the fast 

relaxation of |*Ti '') {i = 2-6) comes from the same spin 
configuration of the initial and final states. As discussed 
above, the relaxation rates of the states with single valley 
state "— " are insensitive to the valley splitting. There- 
fore, the results in the case of small valley splitting are al- 
most the same as the case of large valley splitting. More- 
over, we find the results are also robust against the effec- 
tive well width. 

Similarly, one can calculate the relaxation rates of the 
another set of states with the valley configuration . 
The total relaxation rate of |S'*^"'"^) can also recover the ex- 
perimental data pretty well, where the channel between 
|S'(+)) and |irl+^) is the dominant one. Here, |S'(+)) and 
jiji^H-)^ are the lowest singlet and triplet states of this set 
of valley configuration, separately. As for the set com- 
posed by the states with the valley configuration "m" , 
more triplet basis functions (constructed by two single- 
electron basis functions with the same quantum numbers 
n and I) should be included. This makes the results of 
this set of states different from the other two with sin- 
gle valley index "— " or . However, as the energies 
of the states with single valley state are higher than 
the corresponding ones with valley state "— " , we suppose 
the experimental data by Xiao et al^ corresponding to 
the states with "— " valley index. 
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FIG. 2: (Color online) (a) The lowest four energy levels vs. 
perpendicular magnetic field B± in a single QD. The anti- 
crossing/crossing points are labeled as _B ,, (_), -B„/_%„(-) 

'J- 

and Sg(-)rj,(-)- (b) Transition rates vs. the magnetic field. 
The inset zooms the range near B^^_^^(^) with the rates of 

channels S^'^T^"' and tL"'t^~'. In the calculation, 2a = 
4.344 nm and do = 29 nm. 



B. PERPENDICULAR MAGNETIC FIELD 
DEPENDENCE 



In this part, we turn to the perpendicular magnetic 
field case and choose SiGe/Si/SiGe QDs without loss 
of generality. The relative static dielectric constant is 
AC = 11.9 in this structured We start from the structure 
with 32 monoatomic layers of silicon along the growth 
direction of the quantum well as in the parallel mag- 
netic field case, corresponding to a large valley splitting 
2|AiJ = 0.83 meV. With an electric field 30 kV/cm along 
the growth direction, one obtains the strength of the 
Rashba SOC induced by this electric field oq = =f6-06 m/s 
and that of the IIA term bo = T30.31m/s for the SOC 
elements between the states with identical valley index 

-» 43 



is taken to ensure that the lowest levels are constructed 
only by the basis functions with valley index "— ". 

The first four levels in the QD are plotted in Fig. [2lja) 
as function of the perpendicular magnetic field. They 



-(-)\ 



are labeled as IS'(-)), |T| >) (spin up), \T^ >) and \T 
(spin down), according to their major components. The 
shape of the spectrum can be understood from the 
single-electron spectrum of Eq. (|3]). For example, the 
major component of 15''^^-'), i.e., |5'^'-~^), is composed 
by two electrons in |001— ) state, hence the magnetic- 
field dependence of e^c-) is given by 2Eqq^ approxi- 
mately. Similarly, the magnetic-field dependence of the 



triplet \T^ ^) (|T4 0) can be described by E^ 
{Eqq.^ + i Ez, with Ez representing the Zeeman 

splitting), because this state mainly contains the triplet 

basis \Tq^ ') which involves the single-electron functions 
|001— ) and |0— 11— ). The qualitative analysis still works 
even with the strong Coulomb interaction. It is shown 
that the singlet state |<S''-~^) intersects the three triplet 
levels with the increase of the magnetic field. Since the 
crossing and/or anticrossing points show different proper- 
ties on ST relaxation as discussed above, we now analyze 
the intersecting points. We still denote the two-electron 
angular momentum as L = h + I2, but take the spin 
states (S,Sz) instead by considering the perpendicular 
magnetic field. The SOC Hamiltonian can be rewritten 



-(-) 



001 "1" -^0-11 



P-S-), 
(15) 

with the ladder operations P^ and changing L and 
Sz by one unit, respectively. Here, = i iSy. It 
is clear that a state with {L,Sz) can couple with the 
one with (L ± 1, 5*2 =F 1) due to the Rashba SOC and 
the one with (i ± 1, S'z ± 1) due to the IIA term. Ap- 
proximately, the quantum numbers (L, Sz) of |rl ■*) and 
jS*^"^) are (—1,-1) and (0,0), respectively, according to 
the wave functions of \T^'- ^) and \S^'^-'>). Therefore, the 
IIA term couples these states and an anticrossing event 
occurs at the intersecting point (i3g(_)^(-) ~ 0.85 T), 

where an energy gap pops up (~ 0.2 ^eV). Similarly, the 
Rashba SOC results in the anticrossing between jS''^^) 
and |t|~^) (Bg(_).j,(-) - 2.07T). The intersecting point 

between \S'--^} and Ird"^) {B^^_^^(-) - 1.2 T) is simply 
a crossing point. 

The ST relaxation rates together with the transition 
rates between two triplet states are plotted in Fig.[2l^b), 
which shows that the lifetimes of the excited states are 
extremely long (about four orders of magnitude longer 
than the ST relaxation time in GaAs QD^) and strongly 
depend on the strength of the magnetic field. In the 
vicinities of the crossing and anticrossing points, the tran- 
sition rates show intriguing features. For example, at the 
anticrossing point B^^_^^{-), one finds that all the tran- 



=p"4i Moreover, a large effective diameter do = 29 nm sition rates except the one between |r| ^) and |Tf 
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present either a peak or a valley. According to the pre- 
vious workSjSiiZ' the sharp decrease of the transition rate 

between IS'*-"-') and \T^ '') results from the decrease of 
the emission phonon energy. The origin of the features 
of other channels can be understood from Fig.[3l which 
illustrates the major components of the states around 
Bg(_,^(_,, e.g., (red solid curve), jT^^"^) (blue 

dotted one) and \T^^~'') (green dashed one). One notices 
that when the magnetic field approaches i?g(_)^(-), the 

composition of the \t\_ as well as |Tq ^) (not shown) is 
almost invariant, however, the weig 

in 15''^-') significantly decreases (increases) due to the 
spin-mixing from the SOC. As the component of \T^^ •*) 

in |r| is negligibly small, the weight of jS'^^^^) domi- 
nates the relaxation rate. Therefore, the relaxation rate 
between \T^ and |S'^^-') decreases as shown in the inset 

of Fig.[2jb). However, the composition of |Ti '') varies 
in the opposite way, hence the transition rate between 
Irji and |Ti ') presents a maximum at _Bg(_)^(-) . The 

similar feature of the channel involving |Tq •*) can be in- 
terpreted in the same way. Near the anticrossing point 
_Bg(__)^(-), the physics is quite similar and the transition 

rates of all the channels except one from \T^~^) to |t1"') 
present either a peak or a valley. However, in the vicinity 
of the crossing point _B„, ,_(-) , only the transition rate 

* -'-0 

between jS"*^"^) and [Tg shows a sharp decrease due to 
small phonon energy and other transition rates change 
shghtly, because there is no coupling between \S^~^) and 
iTg"-*) and the components of all the states remain al- 
most unchanged. The variation of the transition rates 
far way from the intersecting points can be understood 
from the dependence of the transition rate on the phonon 
energy as mentioned in the previous subsection. ^'^^ 

To indicate the relative contribution of the LA phonon 
mode to the ST relaxation, we remove the TA mode 
from the calculation, vice versa. The magnetic-field de- 
pendence of the relaxation rate of the channel between 
|t1"^) and IS"'")) is plotted in Fig.Hl One notices that 
the relaxation rate of the TA mode is always larger than 
that of the LA mode. Actually, the calculation of the 
other channels (not shown) also reveals similar conclu- 
sion. The reason lies in the different sound velocities of 
the LA and TA phonons. Since the longitudinal sound 
velocity is about twice as large as the transverse one^ 
the momentum of the LA phonon emission is smaller for a 
fixed phonon energy. As a result, the transition rate due 
to the LA phonon emission process is smaller according 
to Eq. Uni). 

In addition, we investigate the influence of the effective 
diameter do on the ST relaxation. The results are plotted 
in Fig.El One notices that the behavior of the transition 
rates is similar to what obtained above by changing the 
perpendicular magnetic field. Here, an anticrossing point 
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FIG. 3: (Color online) Weights of the major components in 
eigenstates, W [= KColOl' with |0 = l^'"'), |r|-^), It!"'), 
and l^o) = iS^f-') (red solid curve), |r^<"'> ( green dashed 
curve), jr^' ') (blue dotted curve)], vs. the magnetic field 
near the anticrossing point ^3(-)rp(-)- 
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FIG. 4: (Color online) Contributions of the LA and TA 
phonon modes to the relaxation rate between |ri ') and 
IS''"-'). The red dashed and blue dotted curves are for the 
LA and TA modes, respectively. Here, 2a = 4.344 nm and 
do = 29 nm. 

between the singlet and one of the triplets (|t1 ■*)) is also 
observed at do ~ 27.4 nm. In the vicinity of this point, 
we also find the relaxation rate between and |S'(~^) 

is strongly suppressed and the rates of other transition 
channels relevant to these two states show a rapid in- 
crease or decrease too. Therefore, the manipulation of 
the ST relaxation by tuning the dot size is also feasi- 
ble. In the experiment, the dot size can be controlled 
electrically.— 

Finally, we also study the case of small valley split- 
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26 28 
do(iim) 

FIG. 5: (Color online) Transition rates vs. effective dot di- 
ameter do. In the calculation, 2a — 4.344 nm and B± = 1 T. 



ting by taking 27 monoatomic layers along the growth 
direction of the quantum well, where 2|A^^| = 0.35 meV 
according to Eq. ([6]). In this configuration, the SOC 
strengths are unavailable in the literature. We ex- 
tract these parameters according to the results of odd 
monoatomic layers calculated by Nestoklon et al^ and 
obtain oq = T2.28m/s and 6o = ±37.93 m/s for the SOC 
elements between the states with identical valley indices 
"=F", when the same electric field (30kV/cm) as the case 
of large valley splitting is applied. One finds that the low- 
est triplet states (denoted as |t1"^), ITq^"^), and iTj""^)) 
are mainly constructed by the single-electron functions 
|001— ) and |001-|-), in a QD with the effective diameter 
18 nm under a low magnetic field. However, the major 
component of the lowest singlet (|S''~'')) remains in the 
same configuration as the case of large valley splitting. 
Interestingly, the second singlet level (|5*'^™-')), whose 
major component is constructed by the single-electron 
basis functions |001— ) and |001-|-), is almost degenerate 
with |Tq™^), which reveals that the intervalley Coulomb 
exchange interaction is rather smalli^ In this case, no an- 
ticrossing point is observed between the relevant states, 
because of the absence of the SOC element between dif- 
ferent valley states when only the lowest subband is rele- 
vant, as mentioned above. Moreover, we find that the re- 
laxations from the three triplet states to |S'*-~-') are much 
slower than those in the case of large valley splitting be- 
cause these triplets and |5'(~^) are in different sets as 
mentioned above. 



IV. SUMMARY 

In summary, we have investigated the ST relaxation 
in silicon QDs with magnetic fields in either the Voigt 
or the Faraday configuration. Our results in the Voigt 
configuration agree pretty well with the recent experi- 
ment in Si/Si02 QDs. We have identified that the origin 
of the relaxation channel in the experiment is between 
the lowest singlet and triplet in the set with single valley 
eigenstate "— (more likely the "— " one). Besides, we 
also predict the enhancement of the ST relaxation pro- 
cess in the vicinity of the anticrossing point due to the 
SOCs when the magnetic field further increases, which 
should be checked by future experiments. We then fo- 
cus on the ST relaxation in the Faraday configuration 
in SiGe/Si/SiGe QDs and discuss the role of the valley 
splittings. In the case of large valley splitting, the low- 
est levels are all constructed by the eigenstates from the 
lowest valley state. We find that the transition rates 
are about four orders of magnitude smaller than those of 
GaAs QDs due to the weak SOC in silicon. The transi- 
tion rates can be effectively manipulated by tuning the 
magnetic field and dot size. From the magnetic-field and 
dot-size dependence of energy levels, we also observe ST 
crossing/anticrossing points. In the vicinity of the anti- 
crossing point, there exists a small energy gap between 
the singlet and one of the triplet states due to the SOC. 
The transition rates of the channels relevant to these two 
states show a sharp increase or decrease. We show that 
the contribution of the TA phonon mode is larger than 
that of the LA one due to the smaller transverse sound 
velocity. As for the small valley splitting, the eigenstates 
from both valley states contribute. We find the ST re- 
laxation rates in this case are much smaller. 
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Appendix A: G IN COULOMB INTERACTION 

G in Eq. (12) is given by 
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where P";' and Wn'^ come from the lateral and vertical parts of the matrix element (n, /, n^, 7|e''^ ""|rt', Z', n'^, 7'), 
respectively. P is^ 

with C*j ; = ^^^^ (n-i) ^^"^ — * + J + (1^1 + l^'l ~ I'' ~ ^l)/2- sgn {x) represents the sign function and W reads 
WC;7 = K>7|exp(ifc^z)K,7'). 
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